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ABSTRACT 

Context. The Tully-Fisiier relation represents a connection between fundamental galaxy parameters, such as its total mass and the 
mass locked in stars. Therefore, the study of the evolution of this relation in the optical and infrared bands can provide valuable 
information about the evolution of the individual galaxies through the changes found in each band. 

Aims. This work aims to study the Tully-Fisher relation at high redshift in the B, V, R, I, and Ks~bands by comparison with the 
local relations derived from a large sample of galaxies in the redshift range 0. l<z<0.3, processed in the same way, and with the same 
instrumental constraints that the high-redshift sample. 

Methods. Using the large photometric information available in the AEGIS database, we have studied the best procedure to obtain 
reliable k-corrections. Instrumental magnitudes are then k- and extinction corrected and the absolute magnitudes derived, using the 
concordance cosmological model. The rotational velocities have been obtained from the widths of optical lines using DEEP2 spectra. 
At high redshift, this method is found to give better results than using the rotation curve, due to spatial resolution limitations. Finally, 
morphology has been determined via visual classification of the HST images. From the above information, the Tully-Fisher relations 
in B, V, R, I and Ks-bands are derived for the local and high-redshift sample. 

Results. We detect evolution in the B, V and R-bands in the sense that galaxies were brighter in the past for the same rotation velocity. 
The change in luminosity is more noticeable in the bluer bands. This colour evolution, unnoticed in our previous work (Fernandez 
Lorenzo et al. 2009) has been detected thanks to the more reliable k-corrections carried out in this paper, which included photometry 
from B to IRAC bands. The change in the (V-Ks) and (R~I) colours (for a fixed velocity) could be interpreted as an ageing of the 
stellar populations as consequence of the star formation decrease since z=1.25. In addition, we conclude that spiral galaxies could 
have doubled their stellar masses in the last 8.6 Gyr. 

Key words. Galaxies; evolution - Galaxies: fundamental parameters - Galaxies: spiral - Galaxies: kinematics and dynamics 

1 . Introduction This turbulence is important mainly in the most external region 

of the galaxy, where the gravity is lower. However, the optical 
The relation between luminosity and maxim um rotational veloc- emission does not extent so far Then, the contribution of the 
ity of spiral galaxies (JTully & Fish^|l977|) is an important dis- turbulence to the rota t ional veloci ty is negligible when using 



^ ^ tance estimator, and it h as been used for measuring the Hubble these Unes.'Vogt et_alj ([T99l 11997^ modeled rotation curves via 

Q constant Hq (e.g.|Tully & Pierce 2000). The study of this relation optical long sfit observations, and they found a modest lumi- 

(3 at different cosmic times can probe evolution in galaxy prop- nosity evolution in the B-band T FR (AMsd iigh z-local)<-0.4 

,—1 . erties such as its total mass, or the relation between dark and at <z> ^0.5). However. IrTx et al.l (,1997,) and lSimard & Pritcherl 

•• • luminous matter (IZwaan et al.lll995h. Moreover, any consistent (|1998|) found a stronger evolution (AM^ss- 1 .5 at z=0.4), model- 

. '^ ■ model of galaxy formation and evolution should be able to repro- ing the kinematics of disc galaxies in a similar way. Subsequent 

^ duce the evolution of the Tully-Fisher relation (hereafter TFR). works have measured the evolution of th e TFR, again findin g 

5^ Then, its study has important implications in determining funda- diff^ erent and ev e n con tradictory results. IZiegler et al.1 ( |2002|) . 

. C3 , mental cosmological parameters, in the stu dy of structure forma- and iBohmetaH (l2004 . find a relation at high redshift (z=l) 

tion, and in the evolution of disc galaxies (|Navarro & Steinmetz| shallower than that measured in local samples, and evidence 

I2QO0). for a luminosity evolution with look-back time of AMb ~ -1 

Traditionally, the internal kinematics of nearby galaxies is magnitudes at redshift z=l. In theu" study, contrary to previous 

measured using 21cm line widths, but the sensitivity of radio findings, they claim that less massive galaxies were brighter in 

telescopes hampers its use at high redshift. Then, optical lines the past (at fixed rotational velocity), while the most massive 

become important in evolutionary studies of TFR, and some au- would f ollow the local relation. However, other groups, such 

thors have stu died the relation between both velocity indicators. as lBamford et al.1 (12006!) . find luminosity evolution but no slope 

For example, iMathewson et al.l (1 1 9921) compared the projected change. Moreover, these authors find differences in luminosity 

rotation velocity measured from Ha rotation curves, with the evolution ranging from -0.2 to -2 magnitudes. 

velocity measured from integrated HI profiles at the 50% of To cope with contributions to line widths other than that 

the maximun flux level. They obtain a difference of 10 km/s, producedby rotational velocity, Weiner et al. (2006a) measured 

which is attributed to the fact that HI widths measure not only kinematic line widths (crii,), and spatially resolved rotation and 

the rotational velocity but the internal galaxy turbulence as well, dispersion profiles. They combined line-of-sight rotation and 
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and found that this combination Table 1. Depth and average image quality of each measurement. 



dispersion (Sq ^=Q.5Y-^g,+o-l^) 
correlates well with the integrated line width, demostrating that 
(TiD is a robust velocity indicator, that allows constructing scal- 
ing relations with velocity for a population of diverse kinematic 
properties ( dispersion and rotation-dominated galaxies). In an- 
other paper. IWeiner et alTl (l2006bl) used crio to study the evolu- 
tion in the B and J-band TFRs, finding a slope evolution in both 
bands, and a larger intercept evolution in the B-band (-1.5 mag 
versus -0.5 mag in the J-band up to z=1.2). This slope change, 
which implies a large evolution for the m ost massive galaxies , 
is i n the opposite sense than that found by IZiegler et al.l (|2002|) 
and 'Bohmetal.' (2004*). Although the B-band TFR derived by 
IWeiner et al. (2006b) using resolved rotation velocities also 
supports an intercept evolution, the sam ple is too smal l and n oisy 
to measure slope evolution. In addition. iKassin et al.l (|2007|) de- 
mostrated that the scatter in th e stellar-mass TFR is lower when 
using So.5 instead of V,,,/. In F ernandez Lorenzo et al.l (|2009|) 
we analysed the optical line widths to study the evolution of the 
TFR in B, R and I-bands, demostrating that all optical emision 
lines can be used for determining galaxy kinematics. Also, we 
constructed our local TFRs from data derived in a similar way 
than the sample of high redshift galaxies, and found evidence of 
luminosity evolution in all three bands for the greatest redshift 
range of our sample, l.l<z<1.3. 

In contrast to the variety of results found in the optical bands, 
there seemed to be a consensus on the abse nce of evolutio n 
in the K-band TFR (IConsehce et al.l 2005; Fl ores et al.ll2006h . 
However, in a recent work, iPuech et al., (i2008i) do find K-band 
evolution in the sense that galaxies had been fainter in the past, 
a result opposite to what is found in the optical bands (for exam- 
ple, JBamforfe^alj__2006) /The disagreement between this work 
and that of iFlores et all (l2006h (both use 3D spectroscopy for 
deriving velocities) is due to the method for correcting the ro- 
tation velocity and, above all, the local relation used as ref- 
erence. Several works have studied the local TFR in the K- 
band finding diff'eren t results, wi th a slope rangi ng from -6.88 
(lHammeretal.ll2007h to -11.3 (iVerheiienll200lh . Then, to ob- 
tain reliable conclusions on the evolution of the infrared TFR, is 
crucial to find the reas on of these disagre ements, fixing the lo- 
cal relation. Recently, Masters et alj(l2008h presented a universal 
caHbration of the TFR in the 2MASS J, H and K-bands. They 
showed that in all three bands the relation is steeper for later- 
type spirals, and they obtained a slope of -10.017 in the K-band 
correcting all galaxies to Sc type. 

In this era of precision cosmology, the increasingly large, 
deep a nd accurate galaxy surveys within which the DEEP2 
project (iDavis et a l. 2003, 2007) provides the higher resolution 
galaxy spectra currently available, allow going one step forward 
on the study of the TFR: its colour evolution. Models of disc 
galaxy formation predict different colour evolution depending 
whether the evolution is mainly due t o collapse or accretion 
(see, for example, Westera et all l2002h . In our previous work 
(iFemandez Lorenzo et al. 1 12009). we studied for first time the 
(R-I) TFR evolution up to z<1.3, but the dispersion prevented 
us from distinguishing between one galaxy model from another 

In the present work, we have increased the sample to the 
largest redshift range 1 . 1 <z< 1 .4 for all the Extended Groth Strip 
(EGS) with DEEP2 spectra available, with the aim to confirm 
or not the prima facie evide nce for evolution found previously 
dFernandez Lorenzo et al. 120 09). Also, we will extend our study 
to the V & Ks-bands, in addition to the B, R & I-bands, in order 
to study the most p robable galaxy formation model, since the 
IWestera et alj (l2002h models predict larger colour differences in 
these bands. The local sample (0.1<z<0.3) has been extended to 
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the whole EGS as well. Using this new local sample, we expect 
to set a more reliable TFR slope at z=0 to compare with our 
high-redshift sample. 

This paper is organised as follows. In Sect. 2, a description 
of the data and the sample selection criteria are provided. The 
study of k-corrections and derivation of the absolute magnitudes 
and rotation velocities are described in Sect. 3. The results are 
presented in Sect. 4, whereas the last two sections provide the 
discussion of the results and conclusions. Throughout this arti- 
cle, the concordance cosmology with Qao - 0.7, Q^o - 0.3 and 
Ho = 70 km s"' Mpc"' is assumed. All magnitudes are in the 
AB zero-point system. 



2. Data & Sample Selection 

The sample consists of galaxies in the Groth Strip Survey 
(GSS) sky region. The baseline for spectroscopy target pre- 
selection were the galaxies for which DEEP2 spectra (Data 
Release 3, DR3) in this field were available in the re dshift ranges 
0.1<z<0.3 and l.l<z<1.4. The DEEP2 project dDavis et all 
I2OO3L l20()7h is a suryev us ing the DEIMOS multi-object spec- 
trograph jFaber et al.ll2003h in the Keck telescopes, to study the 
distant Universe. The grating used was the 1200 1/mm one, cov- 
ering a spectral range of 6500-9100 A with a dispersion of 0.33 
A/px, equivalent to a resolution R=/l/A/l=4000. Despite that ID 
and 2D spectra are available, in this work we only use the inte- 
grated spectra provided by DR3 (Home 1986), due to the lim- 
ited spatial resolution at high redshift. These spectra were ex- 
tracted along the locus of constant lambda from the 2D spectra, 
using the routine do_extract .pro. Therefore, the ID spectra 
provided by the DEEP2 team are corrected by effects of tilted 
slits. 

T he photometric d ata used here are part of AEGIS sur- 
vey dDavis et al.l l2007h . B, R, and I-band photo metry wer e 
taken with the CFH12K mosaic camera (|Cuillandre et alJl200l1) . 
installed on the 3.6-meter Canada-France-Hawaii Telescope 
(CFHT). These magnitudes are includ ed in the pEEP2 photo- 
metric catalogue (Data Release 1, DRh lCoil et alJ2004 and the 
magnitude errors from sky noise and redshift are available too. 
The data in the V-band (F606W) were taken from HST cata- 
logue. As IR photometry we used the Ks-band, which was taken 
with Palomar WIRC (.Bundv et al.1 I2OO6.) . Finally, for high-z 
galaxies, we used data taken with the I nfrared Array Cam era 
(IRAC), on the Spitzer Space Telescope ('Barmbv et al. 2008) in 
the IRACl (3.6 fim) and IRAC2 (4.5 //m) bands'. In Table 1, we 
have represented the limited magnitude and PSF (Point Spread 
Function) of each measurement. B, R and I-band magnitu des are 
alread y corrected for Galactic reddening based on Schleg el et al.l 
dl998h dust maps. The V-band has been corrected following the 
same work. The Galactic reddening for the Ks and IRAC-bands 
is lower (<0.005 mag) than the error for these magnitudes, so 
they have not been considered. 
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Since in the present work we want to concentrate in the study 
of the TFR in V and Ks-bands, we Hmited the sample to galaxies 
with these photometric bands, further selected by restricting the 
sample to galaxies with emission lines in their spectra, neces- 
sary for obtaining the rotation velocity. In addition, for the high- 
redshift sample, we only selected the galaxies with IRACl and 
IRAC2-bands photometry available, since they roughly corre- 
spond to NIR photometry at rest-frame and can provide a more 
reliable k-correction in the Ks-band. 

The inclination angle (i) must be a selection criteria as 
well, since no correlation between magnitude and rotation 
velocity can be observed for inclinations lower than 25° 
(iFernandez Lorenzo et al. Il2009h . The inclination angle was cal- 
culated from the major to minor axes ratio as found in the 
HST catalogue, while the inclination errors were obtained 
comparing with the in clination derived using Sextractor 
dBertin & Arnout5" 1996) in the combined V+I HST images. In 
this way, we estimated the mean error to be + 2.5° for the local 
sample, and + 6° for the high-redshift sample. Galaxies almost 
edge-on will be more affected for extinction, but in our sam- 
ple all galaxies have inclinations lower than 80°. Therefore, the 
inclination of the final sample ranges between 25° and 80°. 

The second selection criteria was the morphology. Full de- 
tails of the morphology classification are d iscussed in our pre- 
vious work ( Fernandez Lorenzo et al. 11200 9) to which the inter- 
ested reader is referred. In order to select spiral galaxies, we per- 
formed a visual classification of every galaxy using HST images. 
The objects were divided into five groups: elliptical/SO (1%), 
spirals (66%), irregulars (5%), interacting (5%), and unknown 
(23% ). To classify our visually unknown objects, GIM2D ( SimardI 
11998 ) was used. The objects with a Sersic index lower than 2.5 
were considered as spirals. The objects visually classificated as 
unknown that could not be fitted with GIM2D (5% of "unknown" 
objects), were discarded. Finally, after applying all these criteria, 
we are left with 128 galaxies in the local sample, and 1 13 in the 
high redshift (1.1 <z< 1 .4) sample. 

3. Data Analysis 

The luminosity and rotation velocity of disc galaxies are the pa- 
rameters involved in the TFR. At high-redshift, several, not ob- 
vious, corrections are necessary to reliably obtain these parame- 
ters. Moreover, some corrections can change dramatically the re- 
sults faking a real evolution. Apart from the Galactic extinction, 
that affects the observed magnitudes of local and high redshift 
galaxies in the same way, the corrections required to obtain the 
absolute magnitudes in each band are the k-correction and the 
intrinsic extinction. In this section, we will analyse the most con- 
venient way to calculate the absolute magnitudes and kinematics 
of spiral galaxies. 

3.1. Rest-Frame magnitudes 

To calculate the k-correction in the B, V, R, I and Ks-bands, 
we need to know the spectral energy distribution (SED) of the 
galaxy. Since the SED is not generally known, or at least known 
with the required photometric accuracy, it is necessary to use 
an appropiate set of templates for reproducing the SED of each 
galaxy. Additionally, the errors in the photometric information 
used to fit the template and its spectral coverage can strongly 
affect the best-fit template and then the k-corrections. In the 
present work, we use the photometry included in the AEGIS cat- 
alogues. Since different instruments have been used to obtain the 
magnitudes in different wavelengths, the same aperture might 



not ensure a consistent fraction of light in each band, due to the 
PSF, the seeing or the pixel scale, that act spreading more or less 
the object. Then, we choose the photometry as close to the total 
magnitudes as p ossible, to fit the SEDs. For B, R and I-bands, 
ICoil et all (|2004|) measured the total R magnitudes from an aper- 
ture that systematically contains the whole galaxy flux, whereas 
B and I magnitudes are corrected to total ones using the colours 
(B-R)i- and (R-I)i • respectively. For the V-band, we choose the 
MAG_BEST in the catalogue, which provides a nearby to total 
magnitude of the object. For the K^-band, we used MAG_AUTO 
because it is available for a larger number of objects than the 
aperture magnitudes, and we confirmed that there is no signif- 
icant difference between aperture and MAG_AUTO photome- 
try due to the larger errors in this band. Finally, for the IRAC- 
bands, we used the MAG_APER photometry in an aperture of 
1.5" (only for the high redshift sample) because aperture cor- 
rections derived from average mosaic PSFs were applied to the 
aperture magnitudes but not to MAG_AUTO, so we did not use 
MAG_AUTO in this case. We checked that there is a good corre- 
lation between the photometry in 1.5" and 2.14" apertures, and 
then the whole galaxy must be inside the aperture of 1.5". 

Since in this work the photometric bands available do not 
match the rest-frame optical magnitudes, and the Ks-band pho- 
tometry is very noisy, we performed a carefuU and sistematic 
study for establishing the most reliable procedure to do the 
k-correction on our data using various sets of templates and 
methods (see Appendix A). We found that the best result is 
obtained from the nonnegat ive linear combin ation of five tem- 
plates based on the Bruzual & Charlon (l2003h stellar evolution 
synthesis codes obtained by kcorrect. However, when a nois- 
ier band is present in the data, the k-corrections calculated 
by this code are not suitable, since the photometric error of 
the observed quantity propagates to the rest-frame magnitude 
(Appendix A). We have found that in this case, the best method 
for the k-correction is obtained from photometric information 
that roughly match the rest-frame band for which we want to 
calculate the k-correction. Then, the rest-frame magnitude in 
the noisy band is better calculated directly from the best-fit tem- 
plate instead of k-correcting the observed data. This is the case 
for our Ks-band that is roughly matched at rest-frame by the 
IRAC2-band. On the other hand, the k-correction in our opti- 
cal bands can only be determined by interpolation, a procedure 
that requires available information at larger and shorter wave- 
lengths than the band that we want to correct. Finally, we have 
calculated the rest-frame magnitudes used in this work convolv- 
ing the rest-frame best-fit kcorrect template with the filter re- 
sponse in each band. 

3.2. Intrinsic extinction 

The absolute magnitudes were calculated from the luminosity 
distance corresponding to the measured redshifts, using the 
concordance cosmology. Finally, the absolute magnitudes 
were corrected for intrinsic extinction. This correction is 
basically based on inclination. In this work we have adopted the 
ITuUv et aP (Il998h procedure, which is valid for local galaxies. 
Since the dust content of the galaxies may have evolved with the 
redshift, our extinction correction can be under or overestimated, 
so we will analize both posibilities in the discussion section. 
According to this method, the extinction A^ as a function of 
inclination / in the A - band, is defined to be: 



A-" = n log (a/b) 



(1) 
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Fig. 1. Example spectra (not flux-calibrated) of the high-z galaxies in our sample. The x-axis represents the rest-frame wavelengths 
in A. The green line is the fit to the [OII]/l/l3727 A double line obtained with dipso, imposing that both lines had the same width. 



with 



75 = -0.35(15.31 + Mb) 



7r 



-0.24(15.91+ Mr) 



7/ = -0.20(16.61+ Mi) 
7^ = -0.045(18.01 + Mk) 



(2) 
(3) 
(4) 
(5) 



where a/b is the galaxy major-to-minor axis ratio. These 
equations are valid for the magnitudes in the Vega zero-point 
system so our magnitudes were converted into Vega-system to 
do the extinction c orrection. For the V- band extinction, we used 
the Calzetti's law dCalzetti et alj|2000h and we obtained A'^" = 
0.8A^-o. 



3.3. Kinematic line widtlis 

Most local velocities have been historically obtained from ra- 
dio measurements, usually from 21c m line widths at 50 per cent 
of the peak inte nsity (for exa mple, iGiovanelli et al.|[l997l) . As 
already known (lBosmairi98lh . the rotation curves are not per- 
fectly flat at large radii. Then, since the observed HI and Ho- 
gas emission do not span the same radii, the velocities measured 
from both lines can provide diff'erent results. In addition, turbu- 
lent motions broaden the HI profile and play a role in the opti- 
cal v ersus radio velocity width determinations (iMathewson et al.l 
11992 ). Moreover, the existence of three types of rotation curves, 
depending on the relation between the maximum velocity and 
the velocity of t he flat region jV erheiien 2001), can complicate 
the comparison. IMathewson et al.i (il992) compared the rotation 



velocity measured from Ha rotation curves with the velocity 
measured from integrated HI profiles at the 50% level. They ob- 
tain a ratio of 2Vc(Ha)/W^Q(HI) ~ 0.94, which is attributed to 
the contribution of the turbulence to the measured HI velocity, 
thus overestimating the rotational velocity obtained using this 
line, as already mentioned. This contribution is more important 
in the most external region of the galaxy where the gravity is 
lower, while the optical emission does not spread until this re- 
gion. Applying diverse corrections. lRix et alj(ll997h found a fac- 
tor of 0.86+0.04 between IVdHa) and W^o(HI), that must be 
considered when comparing our results with those obtained for 
other local samples using 21cm. 

In the DEEP2 survey, the long axis of every individual slitlet 
is oriented along the major axis of each galaxy, with an accuracy 
enough for deriving rotation velocities. Masks in each field are 
generally oriented at 2 different mask PAs, which are -49° and 
41° in EGS. In a given mask, each slitlet can be oriented up to 
+30° with respect to the long axis of the slitmask. Then, from 
the combination of the masks and slits PAs, the maximum slit 
misalignments (6) would be of 1 1° or 19°. We have checked this 
for a subsample of galaxies, in the ACS images, and we found 
good agreement between the slit and the major axis, with dif- 
ferences always lower than 20°. This slit misalignment results 
in a negligible error in the velocity as already pointed out by 
iGiovaneUi et all d 19971) for PA offsets les s then 15°. Also, in 
a more recent work, I Weiner et al. I (l2006al) concluded that spa- 
cially resolved kinematic (2D) depends strongly on slit position 
angle alignment with galaxy major axis, but that integrated line 
width does not. For example, applying the standard cos"' (6) cor- 
rection, the difference in log(Vinax) for a 5=20° would be ~0.012 
for a galaxy similar to the Milky Way at z=0.2, which is lower 
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Fig. 2. Comparison between local TFRs obtained from DEEP2 data (0.1<z<0.3) in B, R, I and Ks-bands, and IVerheiienI (1200 ih 
local TFRs (dashed lines). The solid line represents the linear fit to the DEEP2 data. The parameters of the fit are shown in Table 2. 



than the error for log(Vniax) in our local sample. Then, we chose 
not to apply the slit misalignment correction. 

The optical line used in this work to derive the rotation ve- 
locity is not the same for each sample of galaxies due to the lim- 
ited DEEP2 spectral range. We have measured Ha line for the 
local sample, and the [Oll] AA3727 A doub le line for the high- 
redshift sample. In Fernande z Lorenzo et a l. (2009), we demon- 
strated that both optical emission lines can be used for determin- 
ing disc-rotation velocities with the aim to compare different 
redshifts samples. To calculate the rotation velocity, we used: 



2V - 

^ * max — 



A^c 



^sin(i)(l +z) 



(6) 



where ; is the inclination angle, /lo the line central wave- 
length at z=0, and A.A the line width at 20% of peak intensity. 
Equation 6 is obtained from Doppler effect, applying the cor- 
rections related to cosmology and inclination. We determined 
AA with the Gaussian fitting routines in the Starlink package 
dipso ( Howarth et"anil996i) . which calculate the FWHM (full- 
width half-maximum). The spectral range and resolution of the 
DEEP2 spectra were planned for resolving the [Oil] doublet, 
and all our spectra fulfill this specification. In Fig 1. we have 
represented four cases with different S/N and line blending. To 
measure the [0\\\AA'il21 A double line, we imposed that both 
lines had the s ame width. Then, we corre cted for the instrumen- 
tal width as in lFernandez Lorenzo et al.1 ([2009). Where the pro- 
file is Gaussian, the line width at 20% of the peak intensity can 



be compared with other velocity width measurements such as 
the FWHM and velocity dispersion, cr, using 



FWHM 

2.35 



W20 
3.62 



(7) 



Since our objects have been classified visually as spiral 
galaxies (see Section 2), we assume that all galaxies are rotat- 
ing. However, we cannot prove this assumption with the existing 
data for the high redshift sample. 

4. Results 

4.1. Local Relations 



In [Fernan dez Lore nzo et al.1 (l2009l) . we derived the local TFRs 
fitting the DEEP2 data points in the redshift range 0. 1 <z<0.3. In 
this work, the local sample have been extended to all Groth field 
with ACS data available, with the aim to set a more reliable TFR 
slope at z=0 to compare with our high-redshift sample. Despite 
the bands that we want to study in this work are V and Ks, we 
have derived the local TFRs in B, R and I-bands as well, in order 
to compare with our previous results. 

We use an error-weighted least-squares fitting to the DEEP2 
data points in the redshift range 0.1<z<0.3, to estimate the 
slope a, and the intercept b of the local TFR in each band. A 
m ore detailed explanation of t he fitt ing procedure can be found 
in [Fernandez Lorenzo et al. I (l2009l) . Again, we have adopted 
\og{2V,„ax) as the dependent variable in the fit. This is the so- 
called inverse TFR, which is l ess sensitive to lu minosity incom- 
pleteness bias (IWimcklll994t ISchechtedlI980h . The results ob- 
tained for each band are given in Table 2. In addition, we show 
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Table 2. Parameters of the Tully-Fisher relations obtai ned by fit t ing D EEP2 data in the redshift range 0.1<z<0.3 used as local 
reference. In addition, we show the values calculated bv iVerheiienI (120011) for B, R, I and Ks-bands. 



This work 



Band 


a 


b 


Ctouil 


A 


B 


A 


B 


B 


0.030 ± 0.045 


-0.117 ±0.002 


0.129 


0.25 ±0.38 


-8.55 ±0.15 


1.24 ±0.82 


-8.5 ± 0.4 


V 


-0.026 ± 0.046 


-0.116 ±0.002 


0.126 


-0.22 ± 0.40 


-8.62 ±0.15 






R 


0.004 ± 0.042 


-0.1 14 ±0.002 


0.121 


0.04 ± 0.37 


-8.77 ±0.15 


1.54 ±0.82 


-8.9 ± 0.4 


I 


0.114 ±0.039 


-0.107 ±0.002 


0.119 


1.07 ±0.36 


-9.34 ±0.17 


3.05 ±0.83 


-9.6 ± 0.4 


Ks 


0.316 ±0.040 


-0.098 ± 0.002 


0.119 


3.23 ±0.41 


-10.20 ±0.21 


5.04 ± 0.92 


-10.5 ±0.4 



the values calculated by |Verheiie"nl (1200 ll) for the local "RC/FD 
sample (without NGC 3992)" case, by fitting Vmax (the same 
sample used for comparison in our previous work). In Fig. 2, we 
represented both relations in B, R, I and Ks-bands. As in our 
previous work, the slop e of the TFR in e ach band is very similar 
to that found locally bv lVerheiienI (1200 ih . even for the Ks-band, 
and consistent within his errors. Moreover, our Ks-band TFR 
slope is consistent with that derived by Masters et al. (2008) af- 
ter coiTecting for incompleteness, morphology and luminosity 
dependence. Unlike the slope, we obtained a z ero-point of the 
relations lower than those of IVerheiienI (1200 ih . However, this 
difference is very similar in all bands, so it could be atributted 
to the tru e rotati on velocity versus observed line width relation. 
iRix et al] (1 19971) modeled simulated observations of disc veloc- 
ity fields with true circular velocity Vc, and found < cr >fWc - 
0.6, where the observed < cr > incorporates an average over in- 
clination. This result implies that 2Vf - 0.92 W2o. In our case, 
a ratio of < cr >/Vc ~ 0.4, which provides 2Vc = 1.38 W20, 
would be necessary to explain the difference found in our local 
TFRs. iRix et al] (Il997h found that factors as line profile asym- 
metries, fiber size or inclination effects work towards reducing 
< cr >/Vc. Our spectra are not measured through fibers, but we 
have other effects such as misalignment of the slit, that also re- 
duces < cr >/Vc. Since we are going to compare our high red- 
shift sample with a sample of local galaxies which are measured 
with the same instrumentation and following the same analysis 
procedures, both samples are affected in a similar way by these 
factors, and then the differences found would be meaningful. 



4.2. High-Redshift sample 

In Fig. 3 we have represented B, V, R, I and Ks-band TFRs, for 
the redshift range l.l<z<1.4, to test the evidence of evolution 
found in our previous work but with a larger sample of galaxies. 
In order to examine the intercept evolution we have set the slope 
to the local one of the TFR in each band. The large scatter of the 
high-z relation is noticeably reduced when we limit the sample 
to galaxies with inclination / > 45° (see Fig. 3), which is the 
lower value used in some works (for example, iNakamura et alj 
I2OO6I) . If we fit separately objects with / < 45°, we found lower 
evolution than for all objects together. For galaxies with / > 45° 
the opposite happens, i.e. we found larger evolution in the zero- 
point. For galaxies with low inclination, an error in the inclina- 
tion angle yields larger velocity eiTors than for more inclined 
galaxies. Nevertheless, the error in the extinction will be more 
noticeable for galaxies with high inclination. If the error in the 
extinction were the reason of the difference between both incli- 
nation ranges, then for the Ks-band there would be no change 
between the evolution found for galaxies with low and high incli- 
nation, since the extinction is negligible in this band. However, 
when both inclination ranges were fitted independently for the 
Ks-band, we found the same difference that for the other bands. 



Table 3. Magnitude evolution found in each band assuming a 
constant slope. 



Band 


AM 


Cam 


B 


-1.125 


0.535 


V 


-0.708 


0.547 


R 


-0.646 


0.531 


I 


-0.298 


0.545 


Ks 


-0.213 


0.610 



Therefore, a velocity error due to an inclination error seems to be 
the reason of the disagreement between both inclination ranges. 
Although in terms of the velocity, the results are more reliable 
for galaxies with / > 45°, the evolution in the zero-point of the 
TFRs do not change significantly if the whole sample is used. 
Therefore, we fitted all galaxies together to increase, into the 
bargain, the statistical significance of the result. The scatter in 
the TFR and its evoluti on are im portant to constraint models of 
galaxy formation. .Masters et alJ (:2008.) observed that the scat- 
ter in the local K-band TFR increases with decreasing rotation 
velocity. In our local sample, the scatter is 1.21 mag in the K- 
band, but it is reduced to 0.9 mag when we select the objects 
with log(2Vmax)>2.2. For l og(2Vmax)-2 . 2, th e scatter in the 
K-band TFR obtained from [Masters et all (12008.) is -0.9 mag. 
However, the scatter in our high-redshift sample (1.23 mag) is 
larger than the scatter of ~0.5 mag obtained in [Masters et al.l 
(2008) for log(2Vmax)=2.5, that is the average velocity of our 
high redshift sample. Then, although the scatter does not evolve 
significantly with r espect to our local relat ions (consistent with 
the result found by IConselice et alJ (l2005h for their samples at 
z<0.7 and z>0.7), it is two to three times higher when compar- 
ing wi th t he local relations of Ver heijen (2001), M asters et al] 
(l2008h . or lHammer et alJ (l2007h . This is probably due to mea- 
surement errors larger than the errors for z=0 galaxies. 

As we have assumed that the slope does not change with 
the redshift, we can calculate the evolution of the magnitude for 
each band, via the difference between the high-redshift and local 
TFRs for a fixed velocity. The results are shown in Table 3. For 
the B-band, we found a change in magnitude of AMb=-1.13 at 
z = 1 . 25 , which is larger than 2a, and co nsistent with our previous 
work ( F ernandez L orenzo et al. 2009), and that found by other 
authors (e.g.. lBohrn et al...2004: .Bamford et al...2006.) . However, 



for R and I-band, we found weaker evidences for evolution, that 
likely can be attributed to the more reliable k-coiTections car- 
ried out in this work. The most interesting result is a lower dif- 
ference in the intercept as the band is redder, being the change 
in the zero-point for I and Ks-bands less than Icr. The differ- 
ence in magnitude found for the Ks-band is compatible with 
no evolution, and it is in agreement with the results found by 
Con selice et al. (2005) and Flores et al. (2006). Nevertheless, 
Pue ch et alj (l2668i) found evolution in the Ks-band TFR, in the 
sense that galaxies have been fainter in the past. They attributed 
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Fig. 3. Tully-Fisher relation in B, V, R, I and Ks-bands for the redshift range 1 . 1 <z< 1 .4. Black diamonds are galaxies with inclina- 
tion 45°<i<80°, and grey points are galaxies with 25°<i<45°. Thin and dashed lines represent the local TFR and its 2cr uncertainty 
in the offset. The thick line is the weighted least-squares fit to all points. The fit was performed by fixing the slope to the local one 
in each case. 



the disagreement with the result of lFlores et al.l (l2006l) to both 
the local relation slope used as reference a nd the more accurat e 
measurement o f the rotation veloc ities of iPuech et al.l (l2008l) . 
However i f the Puech et aP (l2008h data are compa red with the 
Verheijen (2001) local relation (the same used by 'Flor es et aLl 
2006), no evolution is found. Therefore, the local TFR slope is 
playing a fundamental role in the evolution of the Ks-band TFR, 
and it will be studied separately in the next section. 



5. Discussion 

Several authors have studied the infrared TFR both in the K and 
in the Ks-bands. To compare with thes e works, we need to know 
the relation between both bands. From lGrocholski & Saraiedinil 



(l2002h . a Bessell K-band is related with the 2MASS Ks-band as 
K = Ks-(-0.044+0.003)-(0.000+0.005)(J-K). Assuming that 
2MASS and Palomar Ks-bands are equivalent, the difference 
between K and Ks-bands is within the errors in our Ks-band 
absolute magnitudes, and therefore we can compare our results 
with those obtained in other works. 



5. 7. The slope in the local Ks-band TFR 

In order to study the discrepancies found by the authors in 
the slope of the local Ks-band TFR, we use the Millenium 
Simulations database (Springel et al. 2005) carried out by the 
Virgo Consortium. From this simulation, the catalogue that 
we used, contained in the millimil database, is the table 
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Fig. 4. Local Tully-Fisher relation in the Ks-band for the Millenium mock galaxies. The grey points are the objects with Mbuige 
< 0.15 Mstar and M k< -16. The dashed li ne is the relation obtained with our local sample and the three dot-dash line is the local 
relation obtained bv lHammeret al.l (l2007h for a sample of SDSS galaxies with log(Vy7fl,) >2.2. The solid line is the fit to the data 
forMA:< -18 (left) andMA:< -21 (right). 



DeLucia2006a, from which we took the data corresponding to 
maximum rotational velocity (v^ax); Johnson B, V, R, I, and K 
absolute rest-frame magnitudes (dust extinction included); red- 
shift; mass of bulge (Mhiiic,R), and the total ma ss in stars (Mstar) 
dPe Lucia & Blaizotll2007riCroton et al.ll2006h . Using the upper 
value of bulg e-to-total (B/T) mass ratio for late-type galaxies 
from Laurikainen et al.l (i2007) . we select the objects with Mbuige 
< 0.15 Mstar as representative of spiral galaxies. In Fig. 4 we 
represent the data for Mk< -16. K-band magnitudes have been 
converted into the AB system using Mjf(AB)=M/s-(Vega)+1.9 
dHewett et al. ' 2006). The slope of the local TFR obtained using 
simulated galaxies seems to change with the magnitude range 
considered. For instance, when we fit galaxies with y\-K< -18, 
we obtain a slope very close to our sample of local galaxies. 
However, for Mji:< -21 , the slope is closer to the K-band TFR 
found by iHammer et alJ (l2007l) restricting the sam ple to SPSS 
galax ies with log(V/;fl,) >2.2 (the local TFR used bv lPuech et alj 
|2008|) . For the other bands, we found the same tendency for the 
slope as that found in K-band. 

If we assume that the slopes of the TFRs change with the 
magnitude range considered, then our local relations cannot 
be compared with our high redshift sample, which consists of 
galaxies brighter than M^: - -20.5. In order to compare our high 
redshift data with a local relation derived for the same magni- 
tude range, we make use of the Millenium data. With this aim, 
we calculate the difference in the zero-point between our lo- 
cal sample and Millenium data by setting the same slope in the 
same magnitude range as in our lo cal one. As it happens in the 
comparison with IVerheiieril (1200 ih local relations, we obtained 
a similar difference in the intercept of the relation for all bands, 
that can be attributed to the ratio between rotation velocity and 
line width. Nevertheless, in this case, the d ifference in the ze ro- 
point is lower, being the zero-point of the IVerheiienI (1200 Ih re- 
lations greater than the Millenium local zero-points. Then, we 
calculated the Millenium local TFR relations for the brighest 
galaxies and we apply the shift in the zero-point necessary to 
compare with our high redshift sample. Finally, we fit the data 
of the high redshift sample by setting the slope to the new local 
one determined using Millenium data. We found a larger evo- 
lution in luminosity for all bands, but for I and Ks-bands, the 
difference in the zero-point of the high redshift relation with 
respect to the Millenium local one remains inside 2cr (where 
sigma is the uncertainty in the local offset), being ~\cr for the 



Ks-band. Therefore, we found a lower luminosity evolution as 
passing from blue to redder bands, consistent with the results 
obtained in the present work. Then the no evolution of the TFR 
observed in the Ks-band is reinforced. 

5.2. The colour evolution of disc galaxies 

The colours of galaxies provide information on their stellar con- 
tent and, through evolutionary models, on the history of star for- 
mation. Then, the colour evolution versus rotational velocity, al- 
lows us to study the change in the stellar content of a galaxy with 
respect to its total mass. In Fig. 5 we represented the colours 
(V-Ks) and (R-I) versus rotational velocity for local and high- 
redshift galaxies. For (V-Ks) we have a larger dispersion, and 
the evolution of ~0.35 mag found by setting the slope to the 
theoretical local (V-Ks) relation, is within the dispersion (less 
than Icr). However, there is a clear change of ~0.33 mag with 
the redshift, in the (R-I) colour at fixed velocity (>2cr). Then, 
the galaxies are redder today, which could be due to an ageing 
of the stellar population as consequence of a star formation de- 
crease in the last 8.6 Gyr. As in the case of the evolution in R 
and I-band TFRs, we found a different result than that found in 
Fernandez Lorenzo et al. (2009), also due to the more reliable 
k-corrections carried out in this work, where we have used ad- 
ditional photometric information. 

The colour evolution of simulated disc galaxies 
wit h redshift has been investigated in some works. In 
Fernandez Lorenzo et al. I (12009) we compared the results 
with lWestera et al.l (l2002l) . which studied the colour evolution of 
disc galaxies in two models of galaxy formation: the accretion 
and the collapse model. Using these models, the change in 
colours observed here would be better reproduced by the 
accretion model, but the predicted evolution in the (V-K) 
colour would be double than that found here, and the change 
in magnitude is opposite to that predicted by the simulations. 
Moreover, in a subsequent paper, Westera et al. (2007) investi- 
gated the effect of the Inicial Mass Function (IMF) on the colour 
evolution of disc galaxies using two different IMFs: Salpeter 
and Kroupa. The evolution in the (g-K) colour derived from 
their work (for the Kroupa IMF without absortion) would again 
be larger than the result found here, assuming (g-K)«(V-Ks), 
mainly due to the change in the K-band of ~1 mag predicted by 
the simulations. But this change in K would lead, according to 
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Fig. 5. Colour (V-Ks) vs. rotation velocity (up) and colour (R- 
I) vs. rotation velocity (down). Black diamonds represent the 
galaxies in the redshift range l.l<z<1.4, while open squares are 
the galaxies in the redshift range 0. 1 <z<0.3. Difference between 
the local (V-Ks) TFRs (up) and (R-I) TFRs (down) (thin lines) 
and its 2cr uncertainty (dashed lines) are shown. Thick lines are 
the fit to the high redshift data setting the slope to the local rela- 
tion in each case. 



their evolution of the M^/Lk, to a growth of a factor ~6 in the 
stellar mass since z=1.25. However, from the Mt/Lj and M,/Li 
given in the same work, this change in the stellar mass would 
represent an evolution of AMr=0.36 and AMi=0.26 from z=1.25 
to z=0, which would imply that galaxies were fainter in the 
past. This result is opposite to that found here, and that found by 
Ivan Starkenburg et al.' (2006) for the R-ba nd. Then, our results 
cannot be reproduced by the simulations of lWestera et alj (l2002l 
I2007h 



5.3. Comparison with other works 



As we discuss in' Fernandez Lorenzo et al.1 (l2009h . the evolution 
found in the optical bands is more li kely due to a change in lumi- 
nosity than to a change in velocity. Cons elice et alJ (1200 5) and 
iFlores et aP (12006) obtained the same resul t than in this work, 
i.e. no evolution in the K-band TFR, whereas lPuech et all (l2008l) 
found a change of the K-band TFR zero-point, which they at- 
tributed to an average brightening of galaxies since z~0.6 by 
0.66+0.14 mag. They explain this brightening as a growth in 
stellar-mass by a factor ~2.5 since z= 0.6, as evaluated fr om the 
evolution in log(Msteiiar/^K) found bv lDrorv et al.l (l2004l) . They 
claim that the result is consistent with the gaseous 0/H phase 
abundance of z~0.6 emission-line galaxies, which represents 
half that found in present-day spirals (Rodrigues et al. 2008). 

However, assuming a change in log(M ,to ;;„r/Ly) with redshift 
in the way described bv lDrorv et alj (l2004ft . but no evolution in 



the K-band TFR as found in the present work, we obtain that 
the spiral intermediate-mass galaxies have doubled their stel- 
lar masses since z=1.25, which agrees rather well with e xpecta- 
tions from stellar mass density studies (iDrorv et al.ll2004l) . From 



the ma ss-to-light ratios in u, r, i and K-bands of lWestera et al.l 
("2007) (they found a similar result for the K-band than that 
of Drory et al. (2004) in the last 9 Gyrs), but considering our 
growth in stellar mass by a factor 2, the evolution in the opti- 
cal bands would be: AM„=-1.3; AM,=-0.79; AM,=-0.67 (for 
a IMF of Salpeter without absorption), which are consistent 
with the result found in the present work. Assuming at least 
the same evolution in the K- band magnitude at z= 1 .25 as that 
found by iPuech et all ( l2008l) at z=0.6, the stellar mass should 
have grown at least by a fa ctor of 3.5 applyin g the evolution 
in \og(Msieiiar/^K) fouud by iDrorv et al] (|2004|) . However, the 
stellar mass should grow by a factor of 2.6 since z=1.25, if 
the change in the gaseo us O/H abundance wer e attributed to all 
gas transformed in stars dRodrigues et alj|2008l) . Given the short 
lives of the most massive stars responsible for the enrichment of 
the ISM (interstellar medium), and the time required for cool- 
ing the gas to form new stars, the change in O/H in the last 8.6 
Gyr will not be only due to a rise in O, but also to a drop in H. 
Then, likely the change in stellar mass derived from the metal- 
licity evolution is an upper limit. 

The stellar mass can also increase by galaxy merger pro- 
cesses or accretion of small galaxies. In the first case, the galaxy 
disc would be destroyed if the mass of the merger were larger 
than ~15% of the parent galaxy mass. The accretion of small 
galaxies, although a probable contribution in increasing the stel- 
lar mass, it is not enough for explaining a growth by a factor 3.5. 
Hence, our result of a double stellar mass since z=1.25 seems 
more plausible and consistent. Likely the increase in stellar mass 
is due to gas transformed in s tars, consistent with h igher SFR 
density at z~l-2 than today dBouwens et alj 120091) . This gas 
could be prese nt either in the disc, the halo or both places. In 
a recent work, iTacconi et alj ( 1201 Ol) found empirical evidence 
that the larger star formation rates at z~l-2 are the consequence 
of large molecular gas reservoirs (both due to continuous gas in- 
flow from the halos, and minor but not major mergers) and not of 
a larger star formation efficiency than that at z~0. Nevertheless, 
doubling the amount of gas at z=1.25, to account for doubling 
the stellar mass at z=0, would result in hi gher extinct i ons in 
the optical bands than those derived from the lTullv et alj ( |1998[) 
equations used in this work, because they have been obtained us- 
ing local galaxy samples. Were this the case, the evolution found 
in the optical bands would be larger than the evolution observed 
in this work. 



6. Conclusions 

We have studied the evolution of the TFR in B, V, R, I and Ks- 
bands, using two sample of galaxies in the GSS, one in the red- 
shift range 0.1<z<0.3 as representative of local galaxies; and 
the high redshift sample in the redshift range l.l<z<1.4. The 
rotation velocity of all 241 spiral galaxies were measured from 
optical lines widths, using DEEP2 spectra. Morphology was ob- 
tained from HST images, and the absolute magnitudes were de- 
rived from a carefully study of the k-correction made with a 
large set of photometric information. The results of this study 
can be summarized as follows: 

We analized four sets of templates and we concluded that 
the SED of a spiral galaxy is better reproduced by the non- 
neg ative linear combina tion of five templates based on the 
Bru zual & CharloJ (l2003h stellar evolution synthesis codes ob- 
tained by kcorrect. However, when the data include a noisy 
band, the k-corrections calculated by this code are not suitable. 
Also, we find that the most reliable k-correction is obtained from 
information in a observed band that roughly match the rest- 
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frame band. In this case, the rest-frame magnitude in the noisy 
band is better calculated directly from the best-fit template. 

When the observed photometry does not match the rest- 
frame bands, it is possible to calculate the rest-frame magni- 
tudes via an interpolation method, as long as enough informa- 
tion at larger and shorter wavelengths is available. However, an 
extrapolation method is not reliable for galaxies at z~1.25 since 
the results can vary as more than one magnitude mainly due to a 
worse fit of the template. 

New local TFR s were constructed incr easing the sample of 
our previous work (iFernandez Lorenzo et al. 2009) in the red- 
shift range 0.1<z<0.3. We obtained a similar difference in the 
zero-poin t of the TFR for all bands when we compare with the 
IVerheiienI (12001 h local relations, that we attributte to the ratio 
between rotation velocity and line width. 

We derived the high redshift l.l<z<1.4 TFRs fixing the 
slope to the local one in each band, to study the evolution of 
the zero-points. We confir m the evolution in the B-band T FR 
found in our previous work dFemandez Lorenzo et al. 1120091) . in 
the sense that galaxies were brighter in the past, and no evolution 
in the Ks-band TFR. Furthermore, for the other bands we find a 
gradual evolution, with a change in luminosity more noticiable 
as the band is bluer 

We study the slope of the local Ks-band TFR compar- 
ing with Millenium simulated data. We obtained that the slope 
changes depending on the luminosity range used in the fit. 
Considering this effect, we calculated new local Ks-band TFR 
from the Millenium data in the same luminosity range than our 
high redshift sample, confirming again no evolution in the Ks- 
band TFR. 

Assuming t he change in t he ste llar mass-to-light ratio 



M[Lk found bv [ Westera et al.l (l2007h . consistent with that of 
iDrorvet aP (|2004"), but no evolution in the Ks-band, we obtain 
that galaxies would have duplicated their stellar mass in the last 
8.6 Gyr Likely, this growth is mainly due to the gas been trans- 
formed into stars, as supported by the change in the gaseous O/H 
abundance found by Rodrigues et al. (2008). From the mass-to- 
light ratios derived by Westera et al. (2007) for the optical bands 
(IMF of Salpeter), and assuming this change in stellar mass, we 
infer an evolution in the optical luminosity similar to that found 
in the present work. However, the luminosity and colour evolu- 
tion simulated by Westera et al. (2007) is not consistent with that 
found here and would imply a growth of ~6 in the stellar mass 
from z=1.25 to z=0. 

Therefore, the galaxies, despite having doubled their stel- 
lar mass in the last 8.6 Gyr, are nowadays fainter and redder 
in the optical bands, which can be attributed to an ageing of the 
stellar populations as consequence of the star formation density 
decrease observed in the last 8.6 Gyr This is supported by the 
change found in the (V-Ks) and (R-I) colours (for a fixed veloc- 
ity), in the sense that galaxies were bluer in the past. 
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Appendix A: Deriving tlie k-corrections 

In this section we present a systematic study using various sets of 
templates and methods to establish the most reliable procedure 
to do the k-correction on our data. 

For this purpose, we have performed the following tests: 

(i) We compare the results obtained fitting different sets of 
spiral templates to the data, to establish the set of templates that 
better fits our high redshift sample (see section A.l). 

(ii) We estimate how good is the rest-frame photometry re- 
produced by the templates, since our photometry does not match 
the rest-frame optical magnitudes, and the Ks-band errors are 
large. With this aim, we calculate the best-fit template of z= 1.25 
galaxies with aknownSED (i.e.: templates of local galaxies red- 
shifted to z= 1 .25), from the photometric information in the same 
bands and using the same weights that we have in our real sam- 
ple. Then, we obtain the differences between the rest-frame op- 
tical magnitudes calculated from the original rest-frame SED, 
and those from the best-fit template (see section A. 2). 

(iii) We compare the accuracy of the k-corrections derived 
using this interpolation method with those obtained using an 
observed-band that roughly match the rest-frame magnitude 
(see section A. 3). 

A.1. Best set of templates 

To evaluate the templates that better fit each of our galaxies, we 
used the code bpz, which implements the bayesian photometric 
redshift method described in Benitez (2000). We restrict the code 
redshift range to 0.95-1.4, where all the galaxies of our sample 
are located. This code was applied using four different sets of 
spiral galaxy templates: 

(i) The templates of SWIRE, formed by 7 spirals rang- 
ing from early to late types (SO -Sdm), and generated with the 
GRASIL code dSilva et al.lll998h . 

(ii) The templates of Poggianti ('Poggi antil 19971) . that consist 
on models of Sa and Sc spirals computed for various redshifts 
from to 3. 

(iii) The lBruzual & CharlotI (l2003h (hereafter BC03) collec- 
tion of galaxy templates used by iTremonti et al.1 (l2003h in the 
analysis of SDSS (Sloan Digital Sky Survey) galaxy spectra. The 
library includes 39 templates, obtained from 13 models and three 
metallicities Z=0.008, 0.02, and 0.05 (see Table A.l for the de- 
scription of the 13 models). 

(iv) Th e resulting templates ob tained from the routine 
kcorrect (Blanton & Roweis' '2007') which are a nonnegative 
linear combination of five templates based on the BC03 stellar 
evolution synthesis code. 

For each galaxy, we calculate the redshifted best-fit template 
as follows: 






/i = 



(i+zY 



(A.l) 
(A.2) 



where Aq is the wavelength and /i,, is the flux density in 
ergs s"' cm"^ A"' of the template at z=0. Then, we normalize 
to the V-band to obtain the SED of each galaxy. As this tem- 
plate is the best representation of the photometry observed, we 
can know what set of templates reproduces better our data by 
calculating the diff'erence between the observed magnitude niohs 
and the magnitude m,^,,,^ obtained from the redshifted template. 



Table A.l. Description of the templates from BC03 and number 
of galaxies (Ng) in our sample that were fitted to each one by 
bpz. Each model was derived for 3 metallicities: Z - 0.008, 0.02, 
0.05. 



Templates 




Model description 




Age 


Ng 


1-3 




Constant star formation 


6 Gyr 


36 


4-6 




Instantaneus-burst 




5Myr 





7-9 




Instantaneus-burst 




25Myr 


1 


10-12 




Instantaneus-burst 




100 Myr 


21 


13-15 




Instantaneus-burst 




290 Myr 


32 


16-18 




Instantaneus-burst 




640 Myr 


7 


19-21 




Instantaneus-burst 




900 Myr 


5 


22-24 




Instantaneus-burst 




1.4 Gyr 


3 


25-27 




Instantaneus-burst 




2.5 Gyr 





28-30 




Instantaneus-burst 




5 Gyr 


1 


31-33 




Instantaneus-burst 




11 Gyr 





34-36 


Exp 


declining SFR Tsfr = 


5Gyr 


12 Gyr 


8 


37-39 


Exp 


declining SFR tsfr = 


9Gyr 


12 Gyr 


28 



in each band. To calculate m,^,,,^, we project the template into the 
filter Q, using the expresion: 



/_^(^i+i ~ ^^l) ^l f.li TqAi 



niQ - -2.5 log 



^(^,+1 - A,) - fAB (v) r, 



V / 



' QAi 



(A.3) 



where Tq^, is the filter transmission in the Q-band at the 
wavelength Ai, /ab (v) - 3.631 x 10"^*^ ergs s"' cm"^ Hz"' is 
the flux density of the AB standard source, and c is the speed of 
light. 

In Fig. A. 1 we have represented the result for all galaxies in 
each set of templates. For our objects, with a median redshift of 
<z>=1.25, B, V, R and I-bands correspond to UV rest-frame 
wavelenghts, while Ks, IRACl and IRAC2-bands coiTespond 
to NIR rest-frame wavelenghts. For all sets of templates, the 
larger diff'erence is found in the Ks-band due to the photomet- 
ric errors in this magnitude (the fits of bpz and kcorrect are 
weighed by the photometric errors). However, there is a clear 
difference between UV and NIR resulting from the comparison 
with the SWIRE templates (Fig. A. la). For this set, the galaxies 
are bluer than the templates. On the contrary, for the Poggianti's 
templates (Fig. A. lb), the opposite happens, i.e. the galaxies are 
redder than the templates. However, with the set of BC03 (Fig. 
A.lc), the agreement between observed-frame bands and the 
magnitudes obtained from the templates, are better. The com- 
bination of the templates of BC03 obtained through kcorrect 
(Fig. A. Id) provides a similar result, but with smaller disper- 
sions. In this case, the template supplied by the code does not 
need to be normalized to V-band. In addition to the best-fit tem- 
plate, kcorrect gives the k-correction for each band used in 
the fit. This k-correction is the difference between the redshifted 
and rest-frame band, both calculated from the best-fit template. 
Then our rest-frame magnitudes would be obtained by subtract- 
ing this k-corrections to the observed magnitudes. However, the 
photometric errors of the observed quantities propagate to the 
rest-frame magnitudes obtained in this way. This effect spoils 
colour evolution studies if one of the observed bands have sub- 
stantially larger eiTors than the others. This is the case of our 
Ks-band data, so that we used another approach for deriving 
the rest-frame Ks magnitudes. Since the observed IRAC2-band 
roughly match the rest-frame Ks-band, and the fit of the tem- 
plate is quite good in the wavelength range corresponding to 
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4 5 






Fig. A.l. Differences between the B, V, R, I, Ks, IRACl and IRAC2-band magnitudes derived from the best-fit template normalized 
to the observed V-band magnitude (mtemp), and the corresponding observed magnitudes (mobs), for our high-redshift sample, a) 
Templates of SWIRE, b) Templates of Poggianti. c) Templates of BC03. d) Templates provided by kcorrect. 



IRAC2 (see Fig. A. Id), we used directly the magnitudes ob- 
tained in this band from the rest-frame best-fit template in- 
stead of k-correcting the data. For the optical bands, the rest- 
frame magnitudes calculated from the best-fit template or using 
the k-corrections of kcorrect are equivalent in terms of the 
TFR, because the average difference between both methods is 
almost null (see Fig. A. Id). Then, as in Ks-band, the best-fit 
templates obtained through kcorrect will be used, but not the 
k-corrections provided by this software. 



A.2. The reliability of the photometry 

Since our photometry does not match the rest-frame optical 
magnitudes, we need to know how reliably kcorrect can re- 
produce the optical SED of a galaxy at l.l<z<1.4 from the 
photometric information available. With this aim, we have red- 
shifted the templates from BC03 to our median redshift, z= 1 .25, 
and we have calculated the magnitude in B, V, R, 1, Ks, IRACl 
and lRAC2-bands. Then, we have introduced this photometry in 
kcorrect to obtain the best-fit template. We have adopted the 
average magnitude error of our data as the standard deviation 
used by the code in every band. Thus, the bands are weighted as 
the real data. Then, we calculated the difference in magnitude be- 
tween the non-redshifted templates from BC03 and the resulted 
template of kcorrect on the rest-frame B, V, R, I and Ks- 
bands. Such differences are represented in Fig. A. 2 and provide 
a first order estimation of the accuracy of the rest-frame mag- 
nitudes obtained in the previous section. As we have not intro- 
duced photometric information that match the optical rest-frame 



wavelength range, the kcorrect template is providing interpo- 
lated magnitudes for these bands. The general tendency for the 
B-band is that kcorrect provide a fainter rest-frame magni- 
tude that the template from BC03 does, while the Ks-band rest- 
frame magnitudes are brighter However, the average difference 
for V, R, and 1-bands is ~0 when considering the full set of tem- 
plates. For the constant star formation model (first 3 templates) 
and for the two exponential models (last 6 templates), the error in 
the SED reconstructed by kcorrect is ~0.1, whereas the worst 
fit is obtained with a simple stellar population. In general, the 
fit of more recent bursts provides optical rest-frame magnitudes 
brighter than the real magnitudes of the templates from BC03, 
while the optical information of the oldest burst reconstructed 
by kcorrect tends to be fainter. For the simple stellar popula- 
tions that fit some of our galaxies, the errors in some cases are 
~0.4. A recalculation of the k-c orrections in our previous work 
jFernandez Lorenzo et al. 1120091) . adding the Ks-band photom- 
etry to the B, R and I-bands used in this work, shows that our 
previous rest-frame magnitudes were overestimated, especially 
in the 1-band. Moreover, the best-fit template of BC03 provided 
by bpz is different if we use only the information in the B, R and 
1-bands, or if we introduce also the infrared information. We re- 
peated the procedure used in the Fig. A. 2, but introducing only 
the optical bands (B, V, R and 1) in the kcorrect fit, so that 
the rest-frame magnitudes are extrapolated rather than interpo- 
lated. The error in the k-corrected magnitudes is now between 
0.5 and 1.5 mag for most of the templates. Then, the extrapola- 
tion method is not suitable for high redshift galaxies. 
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20 30 
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Fig.A.2. Difference between the magnitudes calculated from 
each rest-frame template of BC03 and from the best-fit tem- 
plate (at z=0) calculated by kcorrect from the photometry of 
the redshifted template at z=1.25. See Table A.l for the descrip- 
tion of the templates from BC03. The grey zones represent the 
templates that were fitted by bpz for some of our galaxies. 
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Fig. A.3. Rest-frame Ks-band magnitude calculated from each 
of the templates from BC03 normalized to an IRAC2 magnitude 
of 21.5 in the redshifts z=l.l, 1.25 and 1.4. The description of 
the templates is the same as in Fig. A. 2. 
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A.3. Comparison of methods 

Now, we are going to study the k-corrections obtained using 
an observed band that roughly match the rest-frame magnitude. 
In our redshift range, IRAC2 gives information about the SED 
region corresponding to rest-frame Ks-band. Hence, if we nor- 
malize the redshifted best-fitting template to the IRAC2-band, 
and we calculate the Ks magnitude from this template at z=0, 
we could have more accuracy in the rest-frame Ks-band mag- 
nitudes. The main difficulty using this procedure is the error in 
the choice of the template. In order to quantify this error, we are 
going to compare the k-corrected Ks-band magnitude obtained 



Fig.A.4. Comparison between the k-corrected magnitude ob- 
tained from kcorrect best-fit template and from the best-fit 
template of BC03. For Ks-band (up), the normalization of the 
template from BC03 was done to the observed IRAC2-band. 
For V-band (down), the template of BC03 was normalized to 
the observed V-band. The solid line shows the least-square fits 
with slope 1 . 



from all templates of BC03 redshifted to z=l.l, 1.25 and 1.4, 
and normalized to the same IRAC2 magnitude. In Fig. A. 3 we 
have represented the results. The larger difference between two 
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Fig. A.5. Comparison between the rest-frame B-band (left) and Ks-band (right) magnitudes calculated from the best-fit template 
obtained using the g, r, i, z, Ks, IRACl and IRAC2-bands aperture magnitudes in 1.5" aperture, and those obtained using the B, V, 
R, I, Ks, IRACl and 1RAC2 total magnitudes. The solid lines are the one to one relations. 



templates is ~0.1 magnitudes at z-lA. We did the same proce- 
dure with the SWIRE templates, obtaining a similar result. Then, 
using photometry that roughly match the rest-frame Ks-band, 
the k-correction is not too dependent of the template, mainly 
due to the smaller difference in the infrar ed SEDs between earl y 
and late-types spirals (see, for example. iMannucci et al.ll200lh . 
However, for optical B, V, R, and 1-bands, this does not hold. 
Some of the galaxies have J-band photometry from Palomar 
data, that match the rest-frame V-band magnitude at z=1.25. 
Following the same procedure that the used with Ks-band (Fig. 
A. 3), but normalizing the templates to mj=22.5, we obtained that 
the maximum difference between the rest-frame V-band of two 
different templates was ~0.3. Since the errors in the J-band are 
larger than the IRAC2 errors, this method is not as good for the 
V-band as for the Ks-band k-correction. 

Finally, we apply both methods to our sample of galaxies in- 
stead of redshifted templates. In Fig. A.4 we have represented 
the rest-frame magnitudes calculated from the best-fit template 
of BC03 normalized to the lRAC2-band, versus the rest-frame 
magnitudes calculated from the best-fit template of kcorrect. 
In addition, we have represented the result in the V-band, but 
the template from BC03 was normalized to the observed-frame 
V-band instead the J-band. For the Ks-band, the difference be- 
tween both methods is negligible, but we found a larger scatter 
between the values of rest-frame V-band. Even so, both values 
are going to provide us the same mean result because the aver- 
aged diff'erences are almost null. 



tal magnitudes in the B, V, R, I, Ks, IRACl and IRAC2-bands, 
which were determined as described in the section 3.1 of the pa- 
per The dispersion obtained in the rest-frame B-band is mainly 
due to the interpolation, the same result that we obtained pre- 
viously (Fig. A. 2). For the rest-frame Ks-band, we obtain very 
good agreement between the magnitudes derived from aperture 
and total photometry, since these k-corrections are basically de- 
pendent to the IRAC2-band, that roughly match the rest-frame 
Ks-band. In this case, the k-correction is not too dependent of 
the template, as we explained in the section A. 3 (and Fig. A. 3). 



A.4. Aperture effects 

All the photometry used in this work is part of the AEGIS 
project and has already been corrected for aperture and PSF 
effects. However, the magnitudes derived at the same aperture 
for all bands are not available. Then, we have chosen to com- 
pare the results obtained in this work with those obtained using 
other optical photometry. With this aim we have used the aper- 
ture photometry in the g, r, i and z-bands, which belong to the 
CFHTLS data and are part of AEGIS too. The aperture used by 
this team is 1 .5", enough to cover the whole galaxies in the high 
redshift sample. Then, we have calculated the best-fit template 
using the aperture photometry in the g, r, i, z, Ks, IRACl and 
IRAC2 bands, with an aperture of 1 .5" (for the noisy Ks-band 
we used mag_AUTO). In Fig A. 5 we have represented the rest- 
frame B and Ks-bands calculated this way, and using the to- 



